Physical properties of ferromagnetic-superconducting coexistent system 
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We study the nuclear relaxation rate 1/Ti of a ferromagnetic-superconducting system from the 
mean field model proposed in Ref.14. This model predicts the existence of a set of gapless excitations 
in the energy spectrum which will affect the properties studied here, such as the density of states and, 
hence, 1 jT\ . The study of the temperature variation of 1/Ti (for T < T c ) shows that the usual Hebel- 
Slichter peak exists, but will be reduced because of the dominant role of the gapless fermions and the 
background magnetic behavior. We have also presented the temperature dependence of ultrasonic 
attenuation and the frequency dependence of electromagnetic absorption within this model. We are 
successful in explaining certain experimental results. 

PACS numbers: 71.27. +a, 75.10.LP, 76.60.Jx, 76.60.-k 



I. INTRODUCTION 

In the BCS theory of superconductivity 1 , the conduc- 
tion electrons in a metal cannot be both ferromagneti- 
cally ordered and superconducting. Superconductors ex- 
pel magnetic field passing through them but strong mag- 
netic fields kill the superconductivity (SC). Even small 
amounts of magnetic impurities are usually enough to 
eliminate SC. Much work has been done both theoreti- 
cally and experimentally to understand this interplay and 
to search for the possibility of coexistence between these 
two ordered states. 

In the conventional theory of SC, ferromagnetism(FM) 
and SC compete with each other but in principle it is pos- 
sible for any metal to become a SC in its non-magnetic 
state at a sufficiently low temperature. Even strongly fer- 
romagnetic iron^ undergoes a superconducting transition 
at low temperature under application of sufficient pres- 
sure to bring it to its nonmagnetic state. But interest is in 
simultaneous existence of both of the ordered states. A. 
Abrikoso^s* studied superconductivity with magnetic im- 
purities using the RKKY interaction in which magnetic 
impurities interact with conduction electrons with the 
magnetization considered as an external parameter inde- 
pendent of the superconducting gap, and showed that the 
normal ferromagnetic state has lower energy than the SC- 
FM state and hence coexistence is energetically unfavor- 
able. Fulde and FerreM studied superconductivity with 
a strong exchange field produced by ferromagnetically 
aligned impurities and found that if the exchange field 
is sufficiently strong compared to the energy gap, a new 
type of depaired superconducting ground state will occur. 
Fay and Appel 5 predicted the possibility of a p-wave su- 
perconducting state in itinerant ferromagnets where the 
pairing is mediated by the exchange of longitudinal spin 
fluctuations. They showed that when the ferromagnetic 
transition is approached from cither the ferromagnetic 
or paramagnetic side, the p-wave transition temperature 
goes through a maximum and then falls to zero. Even if 
superconductivity is interpreted as arising from magnetic 
mediation, it was thought that the SC state will occur 



in the paramagnetic phase. But magnetically mediated 
superconductivity was never observed^. Some theories 
predicted the possibility of s-wave SC and FM coexis- 
tent order in the paramagnetic phase but the ferromag- 
netic fluctuation destroys it near the Curie temperature. 
Consideration of s-wave superconductivity and ferromag- 
netism has been carried out by Suhl 7 and Abrikosov 8 in 
which the ferromagnetism is due to localized spins. There 
had been no experimental observations of coexistence un- 
til it was found in the ferromagnetic metal UGe^. The 
coexistence has also been shown to exist in ZrZn^^ and 
URhGe 11 . Experiments on these three materials show 
that the same electrons are responsible for both ordered 
states. But still the exact symmetry of the paired state 
and the dominant mechanism responsible for the pairing 
is in question. Although most authors believe there is 
triplet superconductivity in these materials, the possibil- 
ity of s-wave superconductivity cannot be denied. 

Blagoev et alMili studied a weak ferromagnetic Fermi- 
liquid and showed that s-wave superconductivity is possi- 
ble and favored on the ferromagnetic side. With similar 
thought, Karchev et al.— developed an itinerant ferro- 
magnetic model in a mean field approach in which the 
magnetic electrons are also the one responsible for the 
formation of the Cooper pairs. Following this model, two 
of the authors^ of this paper studied the specific heat 
and compared it with the experimental data of UGe-i- 
They have shown that this model exhibits the quantita- 
tive behavior of the specific heat. The phase diagram 
shown by them is similar to that found in VGei ex- 
perimentally. So there is enough room to believe the 
existence of s-wave superconductivity in the coexistence 
state. Motivated with these thoughts, in this paper we 
are going to present the results of calculations of the den- 
sity of states, nuclear relaxation rate, ultrasonic attenu- 
ation and electromagnetic absorption in the model given 
by refii^. 

We refer the reader to refii^ for the details, and here 
we look at the mean-field Hamiltonian obtained from a 
model Hamiltonian by the standard mean-field proce- 
dure, 
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The diagonalization of this Hamiltonian using a Bo- 
goliubov transformation yields, 



Hmf = E + J2 { E P a hp + E pP?Pp 
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The quasiparticle energy dispersion relations are, 
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The final step is to minimize the free energy to produce 
the mean-field equations. This results in a set of two 
coupled equations in M and A that will be solved self- 
consistently below. For M we find, 
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and for A, 
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The two order parameters, M and A, have dependen- 
cies such as M = M (g, J, T) and A = A(g, J, T). More- 
over they are coupled with each other through the dis- 
tribution functions which are also functions of M and 
A. From the numerical solutions of these two coupled 
equations (l- r >l(>[) ,we will get the values of M and A to be 
used for the calculation of physical parameters we are 
interested in. We want to clarify that all of the future 
results are derived strictly from the dispersion relations 
and the mean-field equations only, with no other assump- 
tions made about the coupling strength limits or small 
magnetization. In what follows we have studied only the 
case when JM > 2A by which we mean that there al- 
ready exists weak ferromagnetic order in which arises su- 
perconductivity. 
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FIG. 1: Energy spectrum of excited fermions. The upper 
curve corresponds to alpha fermions and lower one is for beta 
fermions. The upper horizontal line is at JM/2 + A and the 
lower horizontal line is at JM/2 — A. The gap is 2A. Beta 
fermion excitations are also associated with the zero energy 
excitations called gapless fermions. 



II. ENERGY SPECTRUM AND DENSITY OF 
STATES 

If we compare the quasi-particle energy spectra 
Eqs.(3,4) with the BCS energy spectrum E = 

<ijQ + |A| 2 , we see that the modification is due to the 
presence of the ferromagnetic energy. The plot of the en- 

2 

ergy spectrum as a function of £ p = ^~ — fi, the energy 
of an excited particle above the Fermi level, is shown in 
Fig.l. It is very clear that the gap is not at the Fermi 
level, rather it is pushed up. Now the energy spectrum 
is symmetric around E = contrary to the BCS case 
where it was symmetric around E = 0. The fermions 
which follow the different energy dispersions have differ- 
ent properties. The fermions which follow the E^ energy 
spectrum arc BCS-like excitations, and will be called 
alpha fermions hereafter, and have a gap in the energy. 
Those fermions which follow the energy spectrum do 
not have a gap and will be called beta fermions or the 
gapless excitations hereafter. The presence of these gap- 
less fermions will change the physical properties of the 
system. The maximum of the energy of beta fermions is 
JM /2 — A and the minimum of the energy of the alpha 
fermions is JM/2 + A and so there will be a gap of 2 A. 
But at T — 0, the beta fermions fill only up to E = 0, so 
the gap of alpha fermions is JM/2 + A. 

Next, we derived the expressions for the density of 
states for the corresponding fermions. We use the usual 
relation, 



N l (E) = JL J dp5{E - El). 
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where i refers to alpha and beta fermions. Using the 
property of the delta function and solving the equation, 
we get, 
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Both of the expressions converge to the density of nor- 
mal fermions at the Fermi level, N(0) = in the 
corresponding limit. And, for M — 0, the expressions 
converge to to the density of s tates of the BCS model, 
N(E)/N(0) = E/^J(E 2 - A 2 ). Both of the expressions 
are the same mathematically, the only difference lies in 
the fact that the energy for gapless fermions ranges from 
—oo to JM/2 — A, including zero obviously and that for 
alpha fermions ranges from JM/2 + A to +oo. The result 
of the calculation of normalized density of states with re- 
spect to density of states of normal state fermions at the 
Fermi level is as shown in Fig. 2. 

As we have mentioned above, both expressions for the 
density of states are identical. If we plot both the expres- 
sion from — oo to +oo, the plot will be same. So we can 
just use any one of the density of states expressions wher- 
ever needed, preserving the limit of energy ranges suit- 
able for the corresponding fermions. In the above plot, 
the left most curve corresponds to the beta fermion den- 
sities and the right most is for the alpha fermions. The 
gap is pushed up (compared to the BCS gap) entirely in 
the positive energy side giving rise to a finite density for 
the gapless fermions at the original Fermi level. The den- 
sity of states is enhanced at JM/2 — A for beta fermions 
and at JM/2 + A for alpha fermions and has no density 
in the gap of 2 A. Looking at the energy spectrum and 
the density of states, we can guess that the physical prop- 
erties of the system will be dominated by normal metal 
like behavior due to the presence of gapless fermions at 
the Fermi energy. 

In order to get a good feeling about the type of gap- 
less fermions we are dealing with, it is better to compare 
this with gapless fermions discussed in the context of a 
superconductor with a magnetic impurity. The origin of 
gapless fermions in FS (our model) and a magnetic im- 
purity superconductor is different. In magnetic impurity 
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FIG. 2: Density of states as a function of Energy of excited 
fermions in our model. There is finite density of gapless 
fermions. The density of state is diverging at JM/2 + A 
and JM/2 - A. The gap is 2A. 



superconductors for zero impurity concentration, there 
will be full fermions condensation, full gap and diverging 
density of states at the edge of gap energy. But an in- 
crease in impurity concentration brings de-pairing effects 
in Cooper pairs. The interaction of fermions with the im- 
purity triggers a flip of the spin of the Cooper pair giving 
rise to what is called a pair breaking effect. The fermions 
resulting from the pair breaking mechanism will start to 
fill in the gap. The diverging density of states starts to 
smear out resulting in a finite density of states for all 
the energy. At certain impurity density the gap will be 
fully populated giving rise to gapless fermions. The zero 
energy gap does not necessarily mean that there is no 
superconductivity because the order parameter will not 
be zero together with the gap for that impurity concen- 
tration. But the transition temperature will be reduced 
indicating a suppressed superconductivity. The super- 
conductivity will be completely suppressed at a certain 
critical impurity density which is bigger than the density 
at which the gap will be zerrA So the gaplessness is due 
to the de-pairing effect of the Cooper pair. But in our 
model of FS we see that the presence of gapless fermions 
is the inherent property of FS. Neither the gap should be 
zero nor should the density of states be smeared out for 
the gapless fermions to exist. There is always Fermi sur- 
face where gapless excitations are present no matter what 
the value of the gap is. Their existence is independent of 
each other. 

Even though the origin of gapless fermions is different 
in these two different cases, their presence has a com- 
mon effect in reducing the transition temperature. In our 
model of FS the presence of gapless fermions is necessary 
for the magnetic property of the system and hence will 
reduce the available fermions to pair up. This reduces 
the transition temperature. Mathematically we can un- 
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derstand how the transition temperature is reduced due 
to the presence of magnetism. Let us compare the equa- 
tion for the gap parameter in the BCS and FS cases. For 
the BCS case, 
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where A is the cutoff value. In the weak coupling limit, 
the equation can be solved to get an expression for the 
gap which looks like, 



A = 2Aexp ( ^To> ) 



(11) 



For our model of FS, referring to Eq.( 14)i4, the equation 
(10) above looks like, 
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and in the weak coupling limit, the gap equation reads 
as, 



A = 2Aexp 



-(^ToT+VC^) 2 -^ 2 ) 
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Analyzing these expressions we clearly see that the pres- 
ence of the magnetization reduces the volume of phase 
space that is available for the Cooper pair. This leads to 
a decrease in the gap at T = and hence T c will also be 
reduced compared to the BCS value for given values of 
other relevant parameters. So the increase in magneti- 
zation pushes the gap away from the Fermi surface and 
reduces its magnitude also. 

Next we study some interesting properties of the coex- 
istent system such as nuclear relaxation, ultrasonic atten- 
uation and electromagnetic absorption following the cal- 
culation done by many authors for BCS model>i§iiLi£ii£. 



III. 



TRANSITION PROBABILITIES 



To study transition probabilities we need to find out 
the expressions for the coherence factors appropriate for 
different cases. For that we use the coefficients of the 
Bogoliubov transformations, namely u and v, 
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Then the expressions for the coherence factors for the 
beta fermions and alpha fermions are, 
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for the scattering process, and 
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for the creation or annihilation process. If we express in 
terms of quasi-particles energy explicitly, these equations 
reduce to 
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for the pair creation or annihilation process. It is clear 
now that F is like a matrix with diagonal components 
for intra band transitions and off diagonal components 
for inter band transitions. The scattering process which 
prefers the upper sign is called case I and the process 
which prefers the lower sign is called case II. In the low 
frequency limit these corresponds to ultrasonic attenua- 
tion and nuclear relaxation respectively. 

The effect of the coherence factors will be discussed 
with respect to the properties studied here. We are now 
ready to have an expression for the net transition rate 
jr, for the coexistent state, between energy levels E and 
E' = E + fiLO , and is expressed below as, 
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\M\ 2 J F(A,E,E')N s (E)N s (E')(f(E) - f{E'))dE 

. (23) 

where M is the magnitude of a one-electron matrix ele- 
ment. Since we are interested in the ratio to the normal 
state scattering rate, we do not need to know more about 
the actual value of M . The limit of integration is from 
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— oo to +00 with an obvious understanding that the lim- 
its of the integral have a gap of 2A, from JM/2 — A to 
JM/2 + A, if we are working with finite temperature but 
the maximum limit of the integral will be only zero if we 
are working at zero temperature since f(E) and f(E') 
both will be zero for E > 0. The coherence factors will 
be chosen appropriate to the band involved in the transi- 
tion. f(E) is the usual probability distribution function. 
We will study the ratio of this transition rate with re- 
spect to that of the normal state fermions at the Fermi 
level (j- = |M| 2 iV 2 (0)M. 

In the limit of Tiuj — > there will be transitions due only 
to scattering and we do not need to consider inter-band 
transitions. This applies to nuclear relaxation and ultra- 
sonic attenuation but for electromagnetic absorption we 
can vary the frequency to higher values so that the inter- 
band transition will also contribute to the scattering. 
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IV. ULTRASONIC ATTENUATION 

The relevant matrix elements for treating the attenu- 
ation of longitudinal sound waves have case I scattering 
coherence factors. So in the limit of hut — > the ratio of 
the transition rate (T„/Ti) can be expressed as, 



FIG. 3: Temperature dependence of low-frequency absorp- 
tion process obeying case I scattering coherence factor. The 
decreasing order of the minima is in the decreasing order of 
energy gap. The dashed curve represents corresponding cal- 
culation for BCS case. 
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We did the numerical calculations of this expression. 
The variation of this transition rate ratio as a function 
of temperature has been presented in Fig. 3. 

There is a significant difference in this graph as com- 
pared to the result based on BCS theory. At T ~ 0, there 
is a finite contribution in the scattering rate whereas in 
the BCS case it is zero. The alpha fermion does not con- 
tribute around this temperature because there will be no 
fermions available to contribute for scattering since all 
fermions will be frozen in the low energy region. Mathe- 
matically it is due to the delta function behavior of the 
differential of the distribution function at energy equal to 
zero so that for integration all the higher energy states 
will be irrelevant. So the finite contribution is from 
the gapless fermions only. For the gapless fermion, at 
E ~ 0, even though the density of states is comparable 
to, but slightly greater than the density of states of nor- 
mal fermions at fermi level, the scattering rate ratio is 
close to but less than one. This effect is due to the co- 
herence factor which will be less than one and dominates 
the effect of the density of states ratio. 

As the temperature is increased the scattering rate is 
decreased. It is definitely due to the effect of the coher- 
ence factors. The logic is the following: With the increase 
in temperature, higher energy states which are close to 



the gap will also be relevant and we know there the effect 
of the coherence factors will be bigger. Even if the ratio 
of the density of states is very high in this energy range, 
the product of the density ratio and differential of the 
distribution function will not be enough to overcome the 
decreasing effect of the coherence factors. If the temper- 
ature is still increased, higher energy states will also be 
excited. The delta function then broadens out and the 
coherence factors increase slowly up to 1 for energy less 
than zero but decrease to zero for energy greater than 
zero. The density of states ratio becomes almost con- 
stant for energy less than zero but enhances for positive 
energy, and the alpha fermion starts to contribute but 
very weakly. For the alpha fermion the coherence factor 
has the most negative effect at the diverging edge of the 
density of states. The future of the scattering rate is de- 
termined by the competition between all of these effects. 
The result is that the rate ratio decreases until the tem- 
perature is close to 0.8T C and starts to increase until it is 
one at T c . From the numerical calculation we see that the 
alpha fermion contribution is responsible for the increase 
in the scattering rate. The curve shows similar behavior 
as in BCS at around T c but as we lowered the temper- 
ature it does not go to zero as in the BCS case since in 
the BCS case no fermion would be available for scatter- 
ing at low temperature. Here, however, the beta fermion 
contribution is comparable to normal metal scattering. 
So if we reduced the temperature from T c we see that 
the initial reduction on the scattering rate is due to the 
alpha fermions which are again BCS- like and later, the 
increase is due to the beta fermions which are gapless. 
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FIG. 4: Temperature dependence of the low-frequency ab- 
sorption process obeying case II scattering coherence factors. 
The increasing peaks correspond to increase in the magnitude 
of the gap. The dotted curve corresponds to the BCS case. 



V. NUCLEAR RELAXATION 
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FIG. 5: The effect of magnetization on temperature depen- 
dence of low-frequency absorption process obeying case II 
scattering coherence factors. The reduction in peak corre- 
sponds to increase in magnetization. The interaction param- 
eters are J = 1.01 fixed and g — 1.082, 1.25, 1.5, 2.0 in order 
of reducing peak, where an increase in g corresponds to an 
increase in M. 



The matrix elements for nuclear-spin relaxation by in- 
teraction with quasi-particles have the case II coherence 
factors, which corresponds to the constructive interfer- 
ence in the relevant low-energy scattering process. The 
scattering coherence factor will have the lower sign. In 
the limit of tko — > the ratio of transitioin rate (T n /T{) 
can be expressed as, 
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The variation of this transition rate ratio as a function 
of temperature has been presented in fig. 4. 

Reasoning similarly as above, the scattering rate is al- 
most equal but greater than one at around zero temper- 
ature. Only the beta fermions contribute at low temper- 
ature since there will be no fermions in the upper band 
which can contribute to scattering. The coherence fac- 
tor has maximum positive effect at — A for the beta 
fermion and ^M- + A for the alpha fermion. The initial 
rise in the ratio of the scattering rate, at a temperature 
close to but less than the critical temperature, is due to 
BCS-like alpha fermions. By around 0.8T C the rate is 
maximum and starts to decrease, due to the freezing out 
of fermions in the beta fermion band. The finite scat- 
tering ratio at lower energy is due to the gapless nature 
of the beta fermions. We solved the coupled equations 
[5,6] for M and A for different values of the interaction 
parameters J and g. We refer the reader to refi^ for 



the details. In fig. 4 the scattering rate corresponds to 
the pairs of J and g for which M is small and remains 
almost constant, and A increases faster. And for those 
values of M and A, we observed that the increase in re- 
laxation rate ratio from around 1.35 to around 1.55 is 
due to the increase in the value of the gap. We see that 
the peak is reduced considerably compared to the BCS 
case and it is due to the presence of the finite magnetiza- 
tion. The increased magnetization reduces the density of 
a-fermions, the magnitude of the coherence factors and 
the probability distribution function. We will show the 
effect of high magnetization explicitly in what follows. 
But here we want to draw attention to the fact that we 
observed the coherence peak which is a signature of an 
s-wave pairing mechanism. 

In Fig. 5, we have presented the effect of magnetization 
on the nuclear relaxation rate. We again solved the cou- 
pled Eqs.(5, 6) for M and A for different values of J and 
g. Fig. 5 corresponds to the scattering rate for pairs of J 
and g for which M is large and increases faster with the 
chosen values of J and g. 

From the graph we can easily infer that even for a fi- 
nite magnetization we can get a finite peak in the relax- 
ation rate. For low magnetization the peak is strong and 
hence we can argue that the superconductivity is domi- 
nant. More fermions are forming pairs than aligning in 
certain direction to give strong magnetization. But with 
the increase in magnetization the peak has been reduced. 
This can be understood by arguing that when the mag- 
netization is increased, the fraction of fermions which are 
giving rise to ferromagnetism will increase with the ul- 
timate effect of decreasing the fraction of fermions that 
form Cooper pairs. We can increase the temperature to 
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get more fermions to form more pairs, but that will still 
hurt the paring since the thermal energy will overcome 
the binding energy of the pairs. Since the coherence peak 
is the inherent property of BCS pairing, once the pairing 
is less effective the peak will be reduced. Nevertheless, 
there are always a finite number of fermions which are 
forming the pairs. So we want to put forward the idea 
that it is natural to have a reduced relaxation due to 
the intrinsic magnetization, but there should in s-wave 
superconductors always survive Cooper pairs to give an 
enhanced relaxation rate, no matter how small it may be. 
In addition, for the time being we believe that the mech- 
anism of superconductivity is s-wave and are optimistic 
about observing the coherence peak in a ferromagnetic 
superconductor, as will be seen in a later section. 
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VI. ELECTROMAGNETIC ABSORPTION 

Unlike the case of nuclear relaxation, it is now possible 
to utilize large enough frequencies to allow the quasi- 
particles to absorb energy, but the absorption process 
will now be different compared to the case of BCS. Now 
there are fermions to absorb any finite energy of electro- 
magnetic wave which was already revealed in the calcula- 
tion of the nuclear relaxation and ultrasonic attenuation 
where the whole curve corresponds to low energy absorp- 
tion. 

The electromagnetic absorption is proportional to the 
real part of the complex conductivity which we can di- 
rectly compute from Eq.(19), with the Fermi function be- 
ing either zero or one at zero temperature. The coherence 
factor will be the one for the pair creation or annihila- 
tion process. We will take the corresponding component 
of coherence factors appropriate for the involved band. 
For different values of incident frequency, the limit of in- 
tegration and the expression to evaluate the absorption 
will be different. For two cases, namely for huj < if — A, 
and if - A < Tiuj < if + A, we use, 
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dE. 



where the limit of integration for the first case is from 
—hu> to 0, where as the limit of integration for the sec- 
ond case is — Jko to — Huj + if — A. These two integrals 
correspond to intra band transition in the beta fermion 
band. In the second case the incident radiation will en- 
counter the effect of the gap. Still, since Tilo < if + A, 
the alpha fermion docs not contribute to absorption. For 



FIG. 6: Frequency dependence of absorption process obeying 
case I coherence factor. 
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where — is the total conductivity and is the sum of con- 
ductivity due to transition from beta to beta 2 ^ — band 
and from beta to alpha band . The limit of integra- 
tion for the first integral is from —Huj to — Tuv + if — A 
and for the second integral is — Jkv+if+A to 0. The up- 
per sign is for case I and the lower sign is for case II. The 
second integral is the contribution of the alpha fermion 
which reduces to a BCS like expression if M goes to zero. 
The first term still gives intra-band transitions and the 
second term gives inter band transitions. It is very easy 
to see that the true gap is only 2A. We evaluated the 
ratio of conductivity for both the cases. The variation 
of this transition rate ratio as a function of temperature 
at fixed value of J,g has been presented in Fig.(6)and 
Fig. (7) for case I and case II respectively. 

Both the cases show almost similar behavior except for 
some peculiarity in each case. For both the cases, the ab- 
sorption rate is never zero because of the presence of gap- 
less fermions. It increases proportionally to the increase 
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FIG. 7: Frequency dependence of absorption process obeying 
case II coherence factor. 



FIG. 8: The nuclear relaxation rate for J = 1.001 and g = 1.1. 
The coherence peak has a 7% enhancement over the normal 
metal value. 



in density of states of the beta fermion if the incident 
energy is enough to excite fermions to a higher energy 
region of the beta fermion band. The case I has lower 
absorption rate than case II if the radiation energy is less 
than the maximum possible value of the beta fermion en- 
ergy because of the coherence factors. Once the incident 
energy is higher than the maximum possible value of the 
beta fermion energy and less than minimum value of the 
alpha fermion energy, the effect of the gap will come into 
play and the absorption rate drops sharply. But this will 
not go to zero since still there is the presence of intra 
band transitions in the beta fermion bands. For energy 
higher than the minimum value of the alpha fermion en- 
ergy the absorption rate will start to increase again and 
reach the value of the normal fermion absorption rate in 
the high energy limit. To be precise, in this energy limit 
the ratio of the rate for case I will be higher and saturate 
out faster than that of case II. If we compare with the 
BCS case, the behavior is similar. The only difference 
in the BCS case is that for case I the ratio has a dis- 
continuous jump when the energy of the electromagnetic 
radiation is bigger than the gap energy and has the ab- 
sorption rate much bigger than that of normal fermions. 
For the case II the ratio starts to increase from zero with 
a faster rate in the BCS case, but here it starts at a finite 
value and rises more slowly. Here, both the inter band 
and intra band transitions will contribute for absorption. 
Another important point, unlike the BCS case, there will 
be no delta function-like behavior of the conductivity cor- 
responding to zero frequency (DC conductivity). In the 
BCS case the lower band will be at an energy which is 
less than zero, hence occupied, and has infinite density of 
states which contribute an infinite conductivity but here 
at zero temperature, the system will never have infinite 
density since an enhanced density corresponds to posi- 
tive energy and will be unoccupied at zero temperature. 
Interestingly enough, both the cases satisfy the optical 



sum rule in the more usual manner even though there 
is a very slight discrepancy in the area bounded by the 
curve, that we guess might be due to the presence of a 
diverging term (density of states) at one extreme limit of 
integration. 



VII. PREDICTION AND DISCUSSION FOR 

UGe 2 

Two of the present authors^ have studied some phys- 
ical properties of U Ge2 using the same models as we 
have used here. In the study of the specific heat of this 
system, they found that their calculation fits well the 
data observed experimentally 20 . The interaction param- 
eters used in that calculation are J = 1.001 and g = 1.1. 
We did the calculation of the relaxation rate ratio with 
respect to the normal state contribution with the same 
interaction parameter. The graph has been presented in 
Fig. 8. It is clear from the figure that the enhancement 
in the relaxation rate ratio is very small, only around 7% 
over the normal metal value. We think it might be hard 
to see this small enhancement experimentally. 

To give a quantitative feeling about the predicted val- 
ues and the experimental values of some properties we 
have presented a table which shows the maximum per- 
centage variation of these physical parameters with re- 
spect to the normal metal values for J = 1.001 and 
.9=1.1. 



property 


predicted 


experimental 


specific heat 


(ACV/CV) = 20% 


(ACV/CV) = 20% 


nuclear relax. 


T„/Ti = 7% 


T n /T x < 15% 


ultras, atten. 


T B /Ti = 1.5% 


T n /T x =? 


electromag. abs. 


<7l/<7 n = 4% 


CTl/o-„ =? 
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FIG. 9: The nuclear relaxation rate for J = 1.001 and g — 1.1. 
The coherence peak is very small. 

It is clear from the data presented in the table that the 
specific heat jump is greatly reduced compared to the 
BCS case with a maximum enhancement of only 20% 
over the normal value. We assume that the maximum 
enhancement in the coherence peak height is less than 
15%, explained below, which is very small compared to 
the more than 100% rise over the normal value in the BCS 
case. The theoretical study of other properties like ul- 
trasonic attenuation and electromagnetic absorption also 
show very feeble variation from the normal state value. 

Experimental data for 1/Ti as a function of T at pres- 
sures around 13Kbar is available in Ref. 21 . To compare 
our prediction with those data, we evaluate our expres- 
sion and present the result in Fig. 9. 

We can see that the coherence peak is very small using 
the interaction parameters that have been used for all of 
the calculations, including the specific heat. In the data 
of Ref.20, the peak is unseen. However, the plot there 
is on a log-log scale, which could reduce its appearance 
based on a reproduction of that data that we have done. 
The point is that the evidence is not clear that there is 
not a coherence peak in the experiments that have been 
carried out so far, and there needs to be a more detailed 
exploration with smaller temperature steps around the 
superconducting transition. We believe there is a very 
small, yet noticeable upturn in the relaxation rate. 



Furthermore, in the experiments by the same group 
on the s-wave superconductor MgB 2 for 1 /Ti— , the plot 
looks incredibly similar to the one for UGe 2 - The peak 
is slightly more noticeable for this well-studied material, 
but there is no background ferromagnetism present to 
suppress it as much as in the case of UGe 2 - 

Another interesting point of our model with respect to 
experiment can be seen in a paper soon to be published 23 
where \jT\T for UGe 2 is studied well below the super- 
conducting transition temperature. It is seen that this 
value holds very constant, independent of pressure and 
temperature. One possible reason surmised for this be- 
havior is the likely presence of low lying gapless exci- 
tations, which is an inherent property of the coexistent 
phase of our model. 



VIII. CONCLUSION 

We closely studied the density of states, ultrasonic 
attenuation, nuclear relaxation and electromagnetic ab- 
sorption for a ferromagnetic superconductor using a 
mean field theoretical approach where superconductiv- 
ity is due to s-wave pairing and the magnetism is due 
to spontaneously broken spin rotation symmetry. Due 
to the finite density of states of gapless fermions, the 
low temperature behavior of the properties we studied is 
dominated by normal metal like behavior. At temper- 
ature less than the superconducting transition tempera- 
ture, our study supports the presence of superconductiv- 
ity of an s-wave nature. Namely, there is the presence of 
a Hebel-Slichter peak. The peak is reduced due to the 
presence of a ferromagnetic background. There is a slight 
decrease in the ultrasonic attenuation rate ratio when the 
temperature is less than the transition temperature. The 
electromagnetic absorption shows not much difference in 
the conductivity for either case. Nonetheless, it shows 
some interesting behavior compared to the result in the 
BCS case. The optical sum rule is followed by both cases 
in a more general way. We are partially successful in 
describing the experimental observation of the tempera- 
ture dependence of the relaxation rate in supercoducting 
UGe 2 . 

We graciously thank to Kotegawa group for providing 
us their experimental data for UGe 2 . This work was done 
with the support of DOE/DEFG0297ER45636. 
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